Reflections from Circular Bends in Rectangular 
Wave Guides — Matrix Theory 

By S. O. RICE 

A method of computing reflections produced by circular bends in rectangular 
wave guides is presented. Tlie procedure employs the theory of matrices. Al- 
though the matrix equations arc (juite simple, a considerable amount of calculation 
is necessary tjeforc quantilalive results may be obtained. Fortunately, the ap- 
proximate formulas pertaining to gentle bends hold surprisingly well for rather 
sharp bends. These formulas are obtained by a limiting'process from the matrix 
equations. The approximate formula for reflection from an H-bend (in which the 
magnetic vector lies in the plane of the bend) generalizes an earlier result due to 
R. E. Marshak. The corresponding formula for the E-bend appears to be new. 

Introduction 

A NUMBER of investigators have studied the propagation of electro- 
-^ *■ magnetic waves in a bent pipe of rectangular cross-section, the bend 
being along an arc of a circle. H. Buchholz , S. Morimoto", and W. J. 
Albersheim have employed Bessel functions to express the field in the bend. 
The form assumed by the field when the radius of curvature of the bend 
becomes large has been obtained by K. Riess and R. E. Marshak who use 
approximations suited to this case. Marshak also obtains expressions for 
various reflection and transmission coefficients. A discussion of the subject 
using rather simple but approximate analysis is given on pages 324-330 of 
a text book'' by S. A. Schelkunoff. The Bessel function approach is also 
sketched in the same section. 

Here we study the disturbance produced when a wave goes around a 
circular bend (of some given angle) in a rectangular wave guide, the guide 
being straight on either side of the bend. Especial attention is paid to the 
dominant mode reflection coefficients gro and d^i corresponding to H-bends 
and E-bends, respectively. As equations (4.2-6) and (4.4-4) show, these 
reflection coefficients (which are of the nature of voltage rather than power 
reflection coefficients) vary inversely as the square of the radius of curva- 
ture of the bend when the bend is gentle. The substance of (4.2-6) has 
been gi\'en by Marshak" for the important case in which only the dominant 
mode is propagated and the angle of the bend not too small. 

When the bend is so sharp that the formulas mentioned above do not 
apply the reflection coefficients maj' be computed from the rather simple 
looking matrix expressions (2.3-3) together with (2.3-4). However, their 
appearance is deceptive and, as is shown by the numerical work in Part V, 
considerable labor is necessary to obtain an answer. 
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The gentle bend formulas were obtained from the matrix equations by the 
limiting process described in Part III. It seems likely that the matrix 
methodj which is similar to the method used in an earlier paper on trans- 
mission line equations, may be applied to other wave guide problems. 
With this thought in mind, the development of Parts II and III has been 
couched in general terms. 

The matrices used in the present theory are of infinite order since the 
guide may support an infinite number of modes of propagation. This fact 
makes it difficult to justify all the steps in our analysis, and we do not at- 
tempt to do so.* Despite this lack of rigor, I believe that the procedures 
given here lead to the correct results since they yield, for gentle bends, 
expressions obtained by Buchholz and Marshak. Moreover, although 
numerical results tabulated in Part V were obtained by using matrices of 
only the second and third order, they indicate a rapid convergence as the 
matrix order is increased. 

PART I 
PROPAGATION OF WAVES IN GUIDE 
1.1 Propagation in a Straight Wave Guide 

Rather general expressions for the electric and magnetic intensities E and 
// in a field are (see pp. 127-128 of Reference ) 

E = —iaifiA -\- -;— grad div A — curl B 

70)6 

(1.1-1) 

H — curl A + T~ grad div B — iunB 

The field is assumed to vary with the time t as e'" , oj is the radian fre- 
quency, n the permeability and e the dielectric constant (for free space 
H = 1.257 X 10'^ henries/meter, t = 8.854 X 10"^^ farads/meter). The 
vector potentials A and B satisfy the wave equations 

V^A = a- A, V'B = <T^B 

V = Laplacian operator (11-2) 

2 2 

In dealing with bends, it is convenient to choose ,1 and B normal to the 
plane of the bend. In our notation, this plane is always taken to be the x, 
2 plane so that A and B are parallel to the y axis. The z axis is parallel 

* Similar questions arise in the rigorous treatment of an infinite set of linear 
equations. A discussion of this subject is given in Chap. TIT of Reference^. 
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to the guide axis and, for the straight guide of the section, the guide walls 
are sections of the planes x- = 0,x= a,y = 0,y = b. 

Thus, a general wave traveling in the positive z direction may be de- 
scribed by the two functions (which represent the magnitudes of A and B) 
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(1.1-3) 



(1.1-4) 



"= 1,2,3, ••• 

where the coefficients g™„ and (it,„ are constants and the plus signs indicate 
propagation in the positive s direction. 

The propagation constant r„.„ is obtained from 

rL = <r'-h (xw/a)- + {irn/b)\ a = i27r/Xo, 
Xo = wavelength in free space. 

Equation (1.1-5) arises when the typical term in (1.1-3) is substituted for A 
in the equation 

^2 * -fit ^^ A 

(1.1-6) 



^^4 a^A ^^4 ^ ^2^ 
a^ 5y dz^ 



This and a similar equation for B are the forms assumed by (1.1-2) for the 
rectangular coordinates of our straight guide. 
The electric and magnetic intensities in the guide are given by 
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(1.1-7) 
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which follow from (1.1-1). 

It is seen that the wave is completely specified by the g^n's and rfmn's. 
These may be arranged as (infinite) column matrices in any convenient 
order. Thus in dealing with (1.1-3) and (1-1-4) we may write 



(1.1-8) 
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In our work we shall consider only those modes corresponding to a fixed 
value of m (or of n) and the order is almost automatically fixed. 

The factors which determine the propagation of the typical terms in the 
summations (1.1-3) and (1.1^) for A and B are 



,+ .-^r, 



a,»n(z) = gmn C 



A.«(2) = ^. e 



(1.1-9) 



The column matrices obtained by arranging these quantities in the same 
order as in (1.1-8) will be denoted by a{z) and &{z). We may write 

«(s) = e-'^" g+, m = r^' d^ 



(1.1-10) 



where exp(—2ro) andexp(— sr^^) are square matrices defined by power series 
each term of which is a square matrix: 



.- = /-£ + 



'J 1 



2! 3! 

/ is the unit matrix and r^, is the diagonal matrix* 



+ 



(1.1-11) 



r„ = 



Tio 















Til 



(1.1-12) 



in which the order of the diagonal elements is the same as the order of the 
elements in the column matrix g+. Similarly F^ is a diagonal matrix whose 
elements are Toi, ro2, Fn, F03, ■ ■ ■ , the order being fixed by d^. When F is 
replaced by Fa in (1.1-11) it is easy to obtain ?„, Fa, etc. and sum the 
resulting series to obtain 



-sFn 



-zFio 











-^Tao 



(1.1-13) 



A similar expression exists for exp(— 3F3). The expression (l.l-lO) for 
a{%) is seen to be true when the square matrix (1.1-13) is multiplied, by 
matrix multiplication, into the column g+. 

It turns out that the field in a circular bend (in a rectangular guide) may 
be represented by a generalization of the foregoing expressions. In this 
generalization, which will be studied in the following sections, the square 
matrices Fo and F^ no longer have the simple form of diagonal matrices. 

* That is, a square matrix in which all of the elements other than those in the principal 
diagonal arc zero. 
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1.2 Propagation in a Circular Bend 

In dealing with a. circular bend we choose cyUndrical coordinates (p, <p, y) 
as shown in Fig. 1. With these coordinates we associate new coordinates, 
shown in Figs. 1 and 2, (.r, y, z) which have approximately the same signifi- 
cance as in the straight guide, s is the distance measured along the axis of 




z=o 



Fig. 1 



!z=o 




Fig. 2 

the guide (defined as the locus of the centers of gravity of the transverse 
cross-sections of the guide), and x and y are the transverse coordinates. 

Let p = pi = (p2 + p3)/2 ^ p3 + "72 be the radius of curvature of the 
guide axis, and let the origin of the polar coordinates be taken at the center 
of curvature. Tiien z is equal to -pi<f> where the minus sign is necessary to 
make (.v, y, s) a right-handed coordinate system. Since the vertical (in 
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Fig. 1) walls are to be specified by x — and .v = a we set .t = p —pi -\- fl/2. 
Thus, the two sets of coordinates are related bv 

p — x -\- pi — a/2 = X -\- Pi 

<P = ~ 2/pi (1.2-1) 

y = y 

where pi, pj and a are constants. 

We again choose A and H in (1.1-1) to be parallel to the y axis. In the 
cylindrical coordinates, 



;we dpdy p d<p p dip iiap, dpdy 



„ I d^A ,dB „ dA , 1 d'/i ,, _ „. 

^v ^ -■ — £-^+7- H^ = - ^- -\- -. — 5-r- (1 .2-2) 

i(atp dipdy dp op lujup oipay 

iii)€ dy' ■iciip, oy^ 

where now, from (1.1-2), A satisfies the wave equation 

liL^_^l+i^_!^ + ^-^ = .^^ (1.2-3) 

P 3pL -^p J p' d<p^ ^/ 

and likewise for B. 
One method of dealing with (1.2-3) which is sometimes used is to assume 

A = c""^ X (sine or cosine function of y) X /(p) (1-2-4) 

where f(p) turns out to be a Bessel function of order p with its argument 
proportional to p. However, we shall proceed in a diflierent direction. 
The change of coordinates (1.2-1) transforms (1.2-2) into 

7W€ dxdy p oz p az m/i axoy 

E, = -ic^A + -^ ^ H,= -iu:eB + i ^^ (1.2-5) 

io}€ By- lojp, ay^ 



mtp dzdy dx dx itaup dzdy 



and (1.2-3) into 



^;4+f^+e?|4+l|l_.M=o. (1.2-6) 

ax' oy^ p'- OZ' p ox 
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where pi is a constant and p ^ x + pi — a/2 is to be considered a function of 
.r. To solve (1 .2- 6) and the corresponding equation for B we assume 

A = ^ a,,,Jz) sin (irmx/a) cos (irny/b) 

(1.2-7) 
m = 1,2, ,^, -■■ ; » = 0, 1,2, ■■- 

/i = S /3„.„(s) cos (irmx/a) sin (irny/b) 

(1.2-8) 
;» = 0,1,2, ••• i » - 1,2,3, ... 

these expressions being suggested by (1.1-3) and (1.1-4). The expressions 
(1.2-5) for the electric intensity show that this choice of .4 and B make its 
tangential component vanish at the walls of the guide. Thus the boundary 
conditions are satisfied. 

In order to determine q!„,„(s) so that the differential equation for A is 
satisfied, we substitute (1.2-7) in (1.2-6). The resulting left hand side of 
(1.2-6) may be regarded as a function, say/(:r, y), of x and y with the a's, 
and their derivatives entering as parameters. We must choose the a's so 
as to make this function zero. Relations which must be satisfied by the a's 
may be obtained by expand ing/(.v, y) in a double Fourier series for which the 
typical term is a coefficient times sin (irmx/a) cos {vny/b), and then setting 
the coefficient of each term to zero. This form of expansion is suggested by 
(1.2-7). However, it should lie mentioned that such an expansion is best 
suited to a function which vanishes at .v = and .r = «, a condition not 
fultilled by/(.v, y) because of the term p~^dA/dx in (1.2-6). This causes no 
real troulile because our region of representation runs only from x = to 
X = a and hence our series is no worse than the Fourier sine series for the 
periodic function (of jieriod 2a) which is — 1 for — o < x < and + 1 for 
< ;v < a. 

To carry out the procedure outlined above, we multiply (1.2-6) (after 
putting in (1.2-7)} by siniirpx/a) cos(irfy/b) and integrate .v from to a 
and y from to b. Using the expression (1.1-5) for r^„ and reducing gives 

-r„fa,A^) + E [F,,..a;j(z) - S,„,a„Az)] = (1.2-9) 

where p may ha\-e any one of the values 1, 2, 3, ■ ■ ■ and the double prime on 
a denotes the second derivative with respect to ::. The P'& and 6"s are 
constants given by 

yj^^^ = (2/a) j {pi/ p') sin iirpx/a) sin (irmx/a) dx, (1.2-10) 

JO 

Su,.. = —2irma~' I sin (irpx/a) cos (wmx/a) dx/p (1.2-11) 
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The evaluation of these integrals is discussed in Appendix I. Thus (1.2-9) 
is the p^^ equation of a set of differential equations to be solved simul- 
taneously for aic{z), oi2C{z), • • • . 

The customary method of solving a set of equations such as (1.2-9) is to 
assume that all the a's vary as e'' so that for each a„u{^) we may write 
e'^gvC- This leads to a set of simultaneous homogeneous linear equations 
for the constants g„,(- In order that these equations may have a solution 
the determinant of the coefficients must vanish. Since the only derivative 
of a„d{z) contained in (1.2-9) is the second, y appears in the determinant 
only as -/. Let 71 , T2 , 73 , • • ■ be the values of 7 which cause the deter- 
minant to vanish and let kij , kij , • ■ ■ be the values of ^i^ , giC , • • ■ cor- 
responding to 7^ = 7j-. The ^'s are determined to within an arbitrary 
multiplying constant which, for the sake of convenience, is chosen so that 
k,j =1. - 

Thus one solution of the differential equation (1.2-6) is 

go 

A - 6"^'' cos (-Kiy/b) ^k^iSiniTrmx/a). (1.2-13) 

m=l 

This particular solution corresponds to the j**" one of the modes (traveling in 
the positive s direction) for which A is proportional to cos {irty/b). 

In much the same way it may be shown that the series (1.2-8) assumed 
for 5 is a solution of equation (1.2-6) (with A replaced by B) provided the 
coefficients j3™„(s) satisfy the set of equations 

-V\i&A^ + i; [Q^^CM - U.^^n^tiz)] = (1.2-14) 

for ^ = 0, 1, 2, ■ ■ ■ and ^ = 1, 2, 3, • - • . Here 

Qpm = (tp/ff) / (p?V) COS {wpx/a) cos {irmx/a) dx (1.2-15) 

JO 

f/j™ = Trmepo"^ / cos (irpx/a) sin {irmx/a) dx/p (1.2-16) 

Jo 

where eo = 1 and e,, = 2 for ^ > 0. These integrals are discussed in 
Appendix I. 

The problem of determinmg the reflection from a bend in a wave guide 
involves considerable manipulation of equations (1.2-9) and (1.2-14). The 
introduction of matrix notation in the manner suggested by the work of 
Section 1.1 for straight guides simplifies this work. Although a„,„{z) and 
?mi>{z) are no longer the simple exponential functions given by (1.1-9), it 
turns out that the column matrices ^(2) and |3(s) are still given by (for a 
wave traveling m the positive 3 direction) by the matrbc expression (1.1-10). 
As mentioned earlier, r„ and T^ are no longer simple diagonal matrices. 
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We now turn to the task of expressing (1.2-9) and (1.2-14) in matrix 
form. 

1 .3 Propagation Constant Matrix for Curved Rectangular Guide 

From this point onward in our investigation of propagation in the rec- 
tangular guide we shall assume A to be proportional to cos (irly/b). Thus 
instead of the general expression (1-2-7) for A we shall deal with the more 
restricted form 



A = cos (irCy/b) ^ a„ai^) sin {irmx/a) 



(1.3-1) 



where ^has one of the values 0, 1, 2, 3, ■ ■ ■ . Since the most general dis- 
turbance may be obtained by the superposition of disturbances of the form 
(1.3-1) no real generality will be lost. 

The introduction of (1.3-1) is suggested by the fact that the set ai^{z), 
a2/(z), - ■ ■ may be determined from (1.2-9) (at least to within arbitrary 
constants of integration) without considering the other a,„n(z)% n ^ I. 
The introduction of (1.3-1) is also suggested by physical reasons. The 
plane of the bend is the z, x plane and there is nothing in the system tending 
to change the field distribution in the y direction. 

Equation (1.2-13) is a special case of (1.3-1). Furthermore the most 
general form of (1.3-1) (corresponding to a wave progressing in the positive 
3 direction) may be obtained by multiplying (1.2-13) by an arbitrary con- 
stant Cj and summing on_/. 

In order to write the set of differential equations (1.2-9) for the a„,?(s)'s in 
matrix form we introduce the infinite matrices 
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v,t 











T^J 











v,c 


Pll 


P12 


-1 


P21 


P22 





i{z) = 



ait{z)' 



(1.3-2) 



S = 



Sn 
Sn 



S22 



where the elements of Tq are obtained by setting n =^ -^ in equation (1.1-5) 
for r„,„, and the elements of P and Q are given by the integrals (1.2-10) and 
(1.2-11). The rules of matrix multipfication then show that (1.2-9) is the 



'■'' element of the matrix equation 

Pa"(s) - (r^ + S)aiz) = 
Premulti plying by F~^ converts this equation into 
a"(z) - rl a{z) = 



(1.3-3) 



(1.3-4) 
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where 

vl = p-'(ro + 5) (1.3-5) 

It may be verified by direct differentiation of the series (1.1-11) defining 
exp {—zV) that a solution of (1.3-4) is 

a{^ = e-'''-g' (1.3-6) 

where, as in (1.1-10) for the straight guide, g+ is a column of constants (of 
integration). However, now r„ is to be obtained by taking the square root 
of the right hand side of (1.3-51, a process which is not easy since it usually 
requires one to obtain the characteristic roots and modal columns of P^ 
(see equation (1.3-10)). 

As far as (1.3-6) being a solution of the differential equation is concerned, 
r„ may be any matrix whose square is given by (1.3-5). We shall restrict 
it as follows: When | — a/p\ becomes small, as in the case of a gentle bend, 
it is seen from (1.2-10) and (1.2-11) that P approaches the unit matrix and 
5. approaches zero. Hence, V a approaches the diagonal matrix To. Ta is 
chosen so that it approaches To, that is, all of the elements in the principal 
diagonal are either positive real or positive imaginary. This makes 
exp(— 2ra) approach the diagonal matrix exp(— sro). With this choice 
of Va the expression (1.3-6) for a{s) corresponds to a wave traveling in the 
positive z direction. 

The various modes of propagation in the bend may be obtained from T „ 
by expressing, in matrix notation, the steps leading to (1.2-13) (which gives 
A for thej"' mode). We assume a(s) to be the column matrix obtained by 
multiplying the column matrix g of constants by the scalar quantity e''. 
Setting this in (1.3-4) gives 

m-Vl)g = (1.3-7) 

where / is the unit matrix. In order that (1.3-7) may have a solution, the 
determinant of the coefficient of g must vanish. This leads to the char- 
acteristic equation* for 7": 

Y^n - Vl\ = (1.3-8) 

The vertical bars denote the determinant of the inclosed matrix. The roots 
71, 72, • • ■ are therefore the latent (or characteristic) roots of r^. K we let 
ifej denote** the column ^obtained when 7 = 7,- in (1.3-7) then 

* See Section 3.6 of Reference^. 

** We choose this notation in order to adhere as closely as possible to that of Refer- 
ence^. Incidentally, the column kj is proportional to the j"' column of k~^ where k is 
the modal row matrix introduced in Section 5.1. 
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and the elements k^, k^, • 
as coefficients in (1.2-13). 



(y)l ~ y\) a, = (1.3-9) 

■ of the modal column kj are the ones appearing 



Equation (1.3-9) and the methods of matrix analysis lead to 

r'„ = ki-iX k~\ v^ = k b], r' (1.3-10) 

where k is the square matrb: whose^ column is kj and [7 ]d, [y]d are diagonal 
matrices having 7,, 7,- as the_; elements in their principal diagonals. The 
representation (1.3-10) certainly holds for the rectangular guide since in 
this case no repeated roots occur. 

In analogy with the expression (1.3-1) for A we shall henceforth deal 
with B in the form 



B = sin {irty/b) 2 &ml{z) cos (Tmx/a) 



(1.3^11) 



where/hasoneof the values 1,2, 3, ■ • • . In much the same way as before it 
may be shown that for a wave traveling along the bend in the positive 
direction the ^^({zys in (1.3-11) are given by 

,3(s) - e"^^ d-^ 



where d^ is a column of arbitrary constants and 



Vl = Q-' {Vl + U) 



(1.3-12) 
(1.3-13) 



In (1.3-12) and (1.3-13) 



Ff, = 



Q = 



Fof 








Tif 





V.J 


(3oo 
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Qio 


Qix ■ 



,3(2) = 






V = 



Uoi Ua2 
f/u U12 

U21 ■ 



^ (1.3-14) 



where the elements of To, Q and U are given by equations (1.1-5), (1.2-15) 
and (1.2-16), respectively. 

1.4 Continuity Conditions at Junclion oJ Straight and Curved Rectangular 
Guides 

Electromagnetic theory requires that E^, Ey, H^ and Hy be continuous in 
crossing a plane 3 = constant which marks the junction of a straight and a 
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curved wave guide (of the same cross-section). Comparison of the first 
equation in (1.1-7) with the first equation in (1.2-5) shows that E^ is con- 
tinuous if (1)^ is continuous and (2) if dB/dz in the straight portion is equal 
(the equality being taken at the junction) to (pi/p) dB/dz in the curved 
portion. Examination of the expressions for the remaining field com- 
ponents shows that all the contmuity conditions are satisfied if, at the 
junction, 

[A in straight portion] = [.,4 in bend] 

and likewise for B. 

Let A in the bend be given by (1.3-1) and let a{z) denote the colxunn 
matrix of coefficients shown m (1.3-2). A in the straight portion may be 
represented in the same way except that a{z) has a simpler form as explained 
in Section 1.1. When these expressions for A are mserted in (1.4-1), both 
sides multiplied by {2/a)sm{irpx/a) after cancelling out the cos (irfy/b), and 
the results integrated with respect to x from to a we obtain relations which 
may be expressed as the matrix equations 

[a{z) in straight portion] = [a{z) m bend] 

where V is the square matrix whose p*' row and jre' column (p, m = 
1, 2, 3;---)is 

Vp,„ = (2pi/a) / sin (rpx/ a) sin iirmx/ a) dx/p, (1.4r-3) 

Jo 

p being equal to pi + a: — a/2. 

By using expression (1.3-11) for B in the continuity conditions, it may be 
shown in much the same way that the column matrix /3(s) given by (1.3-14) 
must satisfy the relations 

\0(z) in straight portion] = [/3(s) in bend] 



["f " '■ " H 



= W — - — in bend 
- l dz J 

where W is the infinite square matrix 

W = T^io T^ii • (1-4-5) 
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whose elemenst are 

Wpm = {(ppi/a) / coslrPx/a) cos (irmx/a) dx/p 

•^ (1.4-6) 

Co = 1, fp = 2 for p > 
Both Vp„, and IFp™ are discussed hi Appendix I. 

PART II 
THEORY FOR A GENERAL WAVE GUIDE 

2.1 Matrix Propagation Constant Jar a Curved Wave Guide of Arbitrary Cross- 
Section 

In Section 1.3 it has been shown that for a curved rectangular wave guide 
there exists a square matrix ?„ (or Tg) which plays the same role in the 
propagation of a wave consisting of many modes as does the propagation 
constant in a simple transmission line. There Ta was obtained from a 
special form of the wave equation which is suited to bends in rectangular 
guides. Here we adopt a different approach with the idea of showing that a 
matrix propagation constant T exists under more general conditions. 

The general theory of wave propagation in tubes shows that a wave 
traveling in the positive z direction may often be represented as 

$ = 12cie-'^^-<Pi{x,y) (2.1-1) 

1=1 

where $ is some quantity associated with the field and is analogous to the 
functions A and B of Part I. In (2.1-1) x and y are transverse coordinates 
and s a longitudinal coordinate, jj is the propagation constant for the 7' 
mode and <pj{x, y) the corresponding eigenfunction. For a circular bend in a 
rectangular wave guide ^j(.v, y) is a combination of trigonometric and Bessel 
functions and 7, is proportional to the order of the Bessel functions. 

We assume that we may find a set of functions 6„,{x, y), m = 1, 2, 3, . . . 
such that every <Pj{Xy y) may be represented as 

v/(-r, y) = E k.„M^, y) (2.1-2) 

m=l 

The usefuhiess of this procedure depends upon our ability to pick a system 
of 6m(x, yYs which is appreciably simpler than the system of <f>i{x, y)'s. In 
the work of Part I dm(x, y) was taken to be the eigenfunction of the typical 
mode of propagation in a straight guide, i.e. the product of a sine and a 
cosine. 
We assume further in (2.1-2) that the square matrix k~^ exists where k„,j is 
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the element in the m*'' row and/'' column of k\ i.e. if a root 7, is repeated, 
say, 5 times there are s linearly independent columns (fey's) corresponding 
to 7/. Substitution of (2.1-2) in (2.1-1) gives 






(2.1-3) 



where 

nM) =flK,iCie-''' (2.1-4) 

Since O^ix, y) is analogous to the product of the trigonometrical terms in 
(1.3-1) or (1.3-11) these equations show that ^(s) plays the same role as 
am^(s) or /9,r.(f(s). Therefore, in accordance with the discussion given at the 
beginning of this section, we wish to show that the column matrbc ju(s) 
(which is similar to a{z) or /3(s)) whose in element is /im(s) may be ex- 
pressed as 

m(2) = e^'V"" (2.1-5) 

In this equation r is a square matrix to be determined and /+ is a column 
matrix of constants similar to g+ or J+. 
The rules of matrix multipHcation and equation (2.1-4) show that 

fx{z) ■= k[e~"']AC (2.1-6) 

in which [exp (— £7)]<j is a diagonal matrix having exp (—£7)) as the / 
element in its principal diagonal and c is the column matrix formed from the 
c,'s. We introduce the column /+ by defining it as ^(0) whence 

/+ =^ kc, c= fe-y+ (2.1-7) 

Incidentally, from (2.1-3), the value of $ at s = is 

*,=o = 'Lfte„.{x,y) (2.1-8) 

171=1 

where /m is the tn'^ element in/*". 
From (2.1-6) and (2.1-7) 

m(s) = k[e-'-']ak-'f (2.1-9) 

In this equation k [exp ( — zy)]d k~^ is a square matrix which may be expressed 
as 

„.„ nl „=.^ «l ^^^^^^^ 

= e ' ., . . 
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Here [7] ^ represents the diagonal matrix havingY/ the;"* term in its principal 
diagonal and 

r = kbUk-' (2.1-11) 

Therefore we have shown that n{z) is of the form (2.1-5) which is what we 
set out to do. 

It is rather difficult to compute T from (2.1-11) using only the above 
definitions of k and 7,- for one must first obtain the functions ifi{x, y). In 
dealing with the rectangular guide it is easier to use equations (1.3-5) and 
(1.3-13) to determine T. 

2.2 Reflection at a Single Junction 

Let a straight wave guide extending from 3= — cotos^Obe joined to a 
curved guide of the same cross-section which extends from s = to s = 00 . 
Let an incident wave 

^i = Jlh„.e-"'"d,Xy,y) ' (2.2-1) 

come in from the left along the straight guide. The kJs are given constants, 
the 6m 's are the modal propagation constants for straight guides (for rec- 
tangular guides they are the r,„n's given by (1.1-5)), and d„,{x, y) is the m* 
eigenfunction for the straight guide (the product of a sine and a cosine for a 
rectangular guide). 

What are the reflected and transmitted waves set up by (2.2-1)? The 
reflected wave is of the form 

^r = T.Ue''-eM,y) (2.2-2) 

m=l 

where the/^'s are to be determined. 
If we assume the representation 

^ = j:nMe„Xx,y) (2.2-3) 

171=1 

to hold for all real values of = then, since 4* = $< + *r for s < 0, equations 
(2.2-1) and (2.2-2) show that 

t^^z) = //„£-''- +/m e'", «; = 1, 2, 3, . . . ; s < (2.2-4) 

Introducing the column matrices t^{z), li, /" and the diagonal matrices 
exp(±sro) where To is a diagonal matrix having 5„ as the m* term in its 
principal diagonal enables us to write (2.2-4) as 

m(s) = e-''" h + e"'"/-, z < (2.2-5) 
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Equation (2.2-5) is more general than the expression (2.1-5) for /i(s) in that 
it contains waves going in both directions, but is more special in that To is a 
diagonal matrix. 

In the curved guide we take m(s) to be given by (2.1-5) , thus 

^(s) = e-'7+, s > (2.2-6) 

where r is a square matrix whose elements are assumed to be known and/+ is 
a column matrix whose elements are to be determined along with those of /~. 
The conditions that the transverse components of the electric and mag- 
netic intensities be continuous at the junction of the two guides are assumed 
to lead to the requirements 

\p.{z) in straight portion] = |ju(s) in curved portion] 

Yd 1 r '^ "I (2-2-7) 

— ^(2) in straight portion = M^ 3- ^(^) ^^ curved portion 

where the quantities within the brackets are evaluated at the junction and V 
is a square matrix whose elements are constants. When the curvature of the 
curved portion becomes small V approaches the unit matrix. For the 
problem at hand (2.2-7) may be written as 

W)U^ = W)],^^ (2.2-8) 

\^^{z)\ = v\^^i^)\ (2.2-9) 

\_dz Jz — \_dz J*-+o 

in which the subscripts 2 = — 0, 3 = +0 refer to the straight and curved 
portions, respectively, of the guide at s = 0. 

The requirements (2.2-7) have been established for the rectangular guide 
in Section 1.4. Their form is also suggested by the conditions that the 
voltage and current be continuous at the junction of two transmission lines. 
Thus if we let }i{z) play the role of the voltage, the current in the first line is 
— ZT^ dn{z)fdz and the current in the second is — Zi" d{i{z)/dz where Zi and 
Z2 denote the distributed series impedances of the two hnes. It is seen that 
this leads to scalar equations which look like (2.2-7), but now V denotes the 
scalar Zi/Za instead of a square matrix. 

Setting the expressions (2.2-5) and (2.2-6) for tj.{z) in the conditions 
(2.2-8) and (2.2-9) gives two equations which may be solved simultaneously 
toobtain/~and/+intermsof /:,ro, rand V: 

(2.2-10) 

-roA + ro/-= -Fr/+ 
/- = (ro+ Fr)-'(ro - VT)h (2.2-11) 

/+ = (Po + VT)-nT(^ (2.2-12) 
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Since f^ and /+ specify the reflected and transmitted waves, respectively, 
they give the answer which we are seeking. 

If the curved guide should extend from 2 = - o) to s - and the straight 
guide from s ^ to s = 00 the response to an incident wave e" h coming in 
along the curved guide would be 

^(s) = e-^^/;+ e""/-, z <0 

(2.2-13) 
m(s) = e-'^^'J^, z >0 

A procedure similar to that used above shows that 

/+ = (Fo + Fr)-i IVIh 
where, instead of condition (2.2-9), we have used 

2.3 Reflection Due to a Bend 

■ Let the guide be straight for - o^ <z < —c and for c < 2 < « , and let 
these two portions be connected by a curved portion in which the longi- 
tudinal coordinate s runs from — c to +c. As in Section 2.2 we take the 
matrk propagation constants for the straight and curved portions to be the 
square matrices To and Y, respectively, and assume an incident wave, 
specified by the column matrix h, to come in from s = — «= . 

The column matrix /x(s) whose in^ element appears as the coefficient of 
&m{x, y) in the representation (2.2-3) for * is now given by 

txiz) = e-'^'h^ e^T, z< -c 

m(s) = (cosh sr)^ + (sinh zV)q, -c < z < c (2.3-1) 

^(.) = .-^"V^ c<z 

In th.£se expressions /~, /+, p, q are column matrices which may be de- 
termined as functions of the known matrices To, T, V and h by substituting 
(2.3-1) in the conditions (2.2-7) which must hold at the junctions z - -c 

and z = c. 

By straightforward algebra similar to that used for the analogous problem 
in transmission line theory we obtain 

-cVo r- -L ^-"^0 f+ ^ r_ r -1- OJVT tanh cT + Tn)"^ Tfil e"^' h 



'f + e-^^»/+ = [-/ + 2iVT tanh cT + To) 'To] e" 
e-'^of- - e-"^^"/"^ = l--^ + 2(l^r coth cT + To)^' To] e'"'" h 



(2.3-2) 
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In these equations the infinite square matrix tanh cV is defined as (sinh cV) 
(cosh cT)~^ and coth cV as its reciprocal, sinh cV and cosh cV may be de- 
fined as power series in cV and may be expressed as combinations of exp (cF) 
and exp (— cr). 

An expression for the column matrix /~ may be obtained by adding the 
equations in {2.3-2). Before doing this it is convenient to introduce the 
two column matrices x and y defined by 

{VcV tanh cT + cFo) x - cVa e^" h 



. (2.3-3) 



- ~ {VcT cot\icV -^ cV^)y = cVfie" h 

where the scalar length c has been introduced to make the various terms 
dimensionless. Each equation in (2.3-3) represents an infinite set of 
simultaneous linear equations to be solved for the elements of x or y. 
Once X and y are known the reflected wave is given by 

/- = e"^" (x + y) - e^'"^' k (2.3-4) 

and the transmitted wave by 

/+ = e''" (x-y) ' (2.3-5) 

PART III 
GENTLE BENDS— GENERAL THEORY ' 

3.1 Limiting Forms Assumed for V and V 

It will be shown in Part IV that for gentle circular bends in rectangular 
wave guides the matrix propagation constant V is such that 

V'' =Vl + F (3.1-1) 



-.2 



where Fo is the square of the matrix propagation constant for the straight 
guide. Fq is a diagonal matrix having 5,„ (which is one of the rLi's, depend- 
ing on the set of modes under consideration, given by (1.1-5)) for the rn^ 
element in its principal diagonal. /^ is a square matrix of infinite order in 
which the elements Fa in the principal diagonal are of order ^^ and the re- 
maining elements Fij, i f^ j are of order ^. Here ^ = a/ pi is the ratio of the 
'guide width to the radius of curvature of the bend. As the bend becomes 
more and more gentle, ^ — » 0. 

The asymptotic expressions given in Appendix I show that, for gentle 
bends in rectangular guides, the square matrbc V which appears in the 
junction conditions (2.2-7) approaches a unit matrix as f — > 0. In par- 
ticular Vii= \ -\- Vii where Va is of order i;^, and V ij, the element in the i^ 
row andy**" column, is of order ^ when i 9^ j. " . 
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Throughout the remainder of Part III we shall assume that T^, F and V 
behave as mentioned above. In addition we assume that there Is no de- 
generacy, i.e. all of the bjs are unequal to each other and to zero. 

3.2 Propagation in a Gentle Bend 

'Here we assume that the elements of To and F in the expression (3.1-1) 
forr- are known. We wish to find the modal propagation constant yy and 
the corresponding eigenfunction tpj(x, y) for they mode. 

After squaring both sides of the coUineatory transformation (2.1-11) 
connecting V and the diagonal matrix {y\d we obtain a relation which may be 
written as k\y'^\d — V^k = 0. The left hand side is a square matrix having 
(7^7 — T^)kj as itsy" column. Here / is the unit matrix and kj is a column 
matrix having kij, ksj, ... as its elements {kj is thej^^ column of k). Thus 
we have a system of simultaneous hnear equations in which the coefficients 
are furnished by the square matrk 7}/ — T- and in which the unknowns are 
kij, kzj-, ■ ■ • . Accordingly, y) is the j^^' latent root of T^ and kj is its cor- 
responding modal column just as for the rectangular guide in Section 1.3, 

In order to apply equations (A2-16) of Appendix II we set >, = y,- and 
M = P so that, from (3.1-1), • 

tiji = 5/ + Fa, Hij = Fii, i^ j (3.2-1) 

Therefore 

y) = S) + Fjj + E' Fj^FsM - b]) (3.2-2) 

k,j = 1, k.j = F.j/{b) - 6l), s ^ j (3.2-3) 

where we have neglected terms of order ^ in (3.2-2) and of order I" in kij, 
s ^ j. The prime on the, summation indicates that the term s = j is to be 
omitted. 

When kij, kij, ... are known the eigenfunction ipi{x, y) may be written 
as a series in 9^{x, y) by means of equation (2.1-2). 

In Section 3.3 we shall need the form assumed by the square matrix V 
tanh cV in a gentle bend. This matrix is used in computing the reflection 
from such a bend, as might be inferred from equation (2.3-3). The formula 
to be used is (A2-18) with u = V-, X; ^ 7^- and with the elements of the 
square matrix k given by (3.2-3). In the diagonal matrix of (A2-18) we set 

/(X;) = X,- tanh cXy = 7,- tanh cy^ = y/tj 

(3.2-4) 
// = tanh cyj 

and for the elements of jfe~' we use (A2-19} together with the line above it. 
When the three matrices on the right of (A2-18) are multiplied out the ele- 
ment (r tanh cr) ij in the i'*' row and^^'' column of T tanh cT is found to be 
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(r tanh cV)ii = (7,(y - JitOka, i 9^ j (3.2-5) 

(r tanh cT)ii = jiti + 2' (jiti - ymtm)ki^k^i (3.2-6) 

where in (3.2-5) and (3.2-6) terms of 0{e) ("order of") and 0(£'), re- 
spectively, have been neglected. This is in line with the fact that the terms 
in the principal diagonals of our matrices must be accurate to within O(f^) 
while the remaining terms need be accurate only to within terms of 0(0). 
The summation with respect to m runs from m = 1 to «> with 
tn = i omitted. 

3.3 ReJlecHon from a Gentle Bend 

When the bend is gentle so that V and T behave according to the descrip- 
tion given in Section 3.1, the matrix expressions for the reflection coefficients 
given in Section 2.3 may be evaluated. The results stated in Appendix II 
for "ahnost diagonal" matrices furnish the principal tools for this work. 

It is assumed that the incident wave coming in along the straight guide 
from the left is $i = exp(~z8p) dp(x, y) and hence contanis only the p^^ 
mode. Comparing this with the general expression (2.2-1) for $i shows 
that hp= l,h^ = 0,m7^ p, and all the elements of the column matrix h are 
zero except the p^^ which is unity. 

We start by writing the first of equations (2.3-3) as 

(r tanh cT + F->ro)flr = K"^ To e"''' h (3.3-1) 

Smce V approaches a unit matrix as $ ^ 0, the element ( K-^ a in the i^^ row 
andy"" column of K~^ is 

{V-')ii = 1 - vii + Z' Vi„ V^i 

m 

where Vii= \ + ''^ihhj ^ 1, 2,3, . . . and the summation with respect torn 
runs from 1 to « with the term iorm = i omitted (as indicated by the prime 
on S). In omitting this term we are neglecting va because it is of order 
t. These results follow from equation (A2-2) of Appendix II. As usual, 
the elements in the principal diagonal are accurate to within 0(^^) and the 
remaining elements to within O(^). 
It follows that V^^ To e'"^" /; = tj is a column matrix whose i*-^ element is 

T/p = (1 - Vpp + 2' Vp„, V„,p) Spc" ", i = p 
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Likewise, the element in the i"' row and y** column of the square matrix 
V~^ To is — VijSj when / ^ j and is 

(1 - Vii + Z' Vi^ V^i) Si (3.3-4) 

m 

when i = j. 

By combining the approximate expressions for the elements of T tanh cV 
(given by (3.2-5) and (3.2-6)) and V-^ To we find that if tu,- denotes the i'-^ 
row and/'' column of T tanh cT + F~* To then 

Mij = Dii ~ Vij8j , i 9^ j 

Uii = yih + E' D,niki,n + 5.(1 - r., + J^' C.i) (3.3-5) 

»i (II 

= 0-. - 5,1',, + E' {D^iki,,, + 5iC„i) 

til 

In these equations we have set 

Oi = 5i + 7,(,-, C„,,- = FimFnii 
Dmi = (Ti^- ^ 7"i'".)^mi = {yiU — ymt„,)Fmi/(d] — 5l) 



(S.S'-d) 



where 7i and ^„,i are given by (3.2-2) and (3.2-3). 

We are now in a position to identify the matrix equation (3.3-1) for x 
with the set of equations (A2-20). The quantity t/,- which appears on the 
right hand side of the i^' equation in (A2-20) is given by (3.3-3). The 
coefficients which appear on the left hand side are thaw's defined by (3.3-5). 
Therefore, from (A2-21), when i 7^ p, 

Xi = -^— [— T',p(ff„ — 5p) — Dip] 



where we have neglected higher order terms and in so doing have replaced 
7j' by the simpler 8j . 

When i = p (A2-21) yields 

Xp = «pp[l?p + 2^ (ll,„.,n-„,p1Jp — tlp„,Upp1J„,)/(Umm''lpp)] (3.3-8) 

In order to combine the second order terms in l/upp with those in the rest 
of the expression for Xp we assume that ap is the major portion of Mpp . 
Then, approximately, 

l/«pp = ffp [1 + SpVpp/a,, — 2J (D^jfkpm + 8pCmp)/<Tp] (3.3-9) 
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This assumption, which is equivalent to assuming that tanh cy^ differs 
appreciably from —1, does not appear to restrict our results since tanh 
C7p is either purely imaginary or else real and positive. 

Substituting the appropriate values in (3.3-8), neglecting higher order 
terms, and using the definition (3.3-6) for Cmp leads to 

Xp = ff~^5„e'^''[l - (1 - 8p/<Tp)vj,j, + XI' (o-pO-m)"'{C„p((rp - 5p)((rm - 5^) 

m 

+ D^piDp^ - cn^kj,^) -\- ((Tp - 5p)I^pmF^p - 5r„Z?„pFp„}] (3.3-10) 

A reduction similar to that used in going from the first to the second line 
of (3.3-7) gives our final expression for Xp 

fip + ypip (1 + Q L 

+ Yl' J. /i"".!™, \/k2 72; {Vpn<5„,F„,p/tp (3.3-U) 

V 5„dp{l ■+■ fm)t5m — Op) 

... - Vr„pSpF^ + Fp„F„.,(5p + S^/tp)(5l - 4)-'} J 

The above expressions for Xi and Xp have been derived from the first of 
equations (2.3-3). The second of equations (2.3-3) determines the column 
matrix y in the same way that the first equation determines x except that 
coth cT now replaces tanh cT. Therefore, we may obtain expressions for 
the elements of y by replacing the i's (where ^ = tanh cyd by their recipro- 
cals in the expressions for the corresponding x's (i.e. m (3.3-7) and (3.3-11). 
The values obtained in this way lead to, when i ^ p, 

/7 = e''{xi-\- yd 

= 57ie'<»(+»p--.-V sinh ^(7^ + yp)[-Vip5„ - Fip/{di + 5p)] 

(3.3—12) 
/I = e''^(x, - yd 

= 57'/»'+*i'-^'-V sinh dyi - yp)[Vip6, - Fip/iSi - 5,)] 

where we have used the expressions (2.3-4) and (2.3-5) for/" and/+. 
When i = p, 

f- = e'Hxp + yp) - €''"'• = -e'-^^-^^'Ui (sinh 2<:7p)/2 + A2] (3.3-13) 

where 

Ai = 2Vpp + (Tp " Op)5p — 2^ Vpm Vi„p -\- ga _ ^2 

„, (cosh 2C7p - g"'""") r-^ ., -2 ,,. „ x^ ^p F ] 

At = 2^ ~ „. , /t2 l2\ [ypmr,„pdn + \'mprp„dp -^ T pm T mpl 

m ^OmOp{0m — Op) 
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The expression for /p may be obtained in the same way but it is sHghtly 
more complicated. 

/t = e\x, - y^) = e^=''--V[i - A, -\- A^ (sinh 2cy^)/2 + A,] (3.3-14) 

where 

A,= (1 - e-'^^0(7p - 5p}V(44) 

' Jim » mp ' mp " pm ^" -Dm " n 



e - T^pm Fmp + ;2 72 + 

m L 



hi - 5^ " (5l - 5'J^ 



. _ ^, (e-'=Tmcosh2co'p - 1) 

m 26m5p(5m — 5p) 



[ 






There are several points we should mention about these formulas for/^ and 
/pi The summations with respect to m run from 1 to m with the term 
m = ^omitted. 7/ and 5/ are the propagation constants of the^th mode 
in the bend and in the straight portion, respectively. The difference 
7p — 5p may be expressed in terms of the /^'s by equation (3.2-2). In the 
course of obtaining (3.3-13) and (3.3-14J relations of the following sort 
were used. 

^p(l + /p)"' = e"''^^(sinh 2c7p)/2 
{U + t\){\ + /p)"'(l + UT' = e~-''-{m^\i 2cyp - e"'^'"') , . 
The term ^3 arises when we subtract (1 — tp) (1 + t^~'' from 
5p/(Sp -f- 7p/p) - 5p/p/(7p + hptp) 

Since 7p — 5p is 0(^^) for a circular bend in a rectangular guide (7^ — 5p) 
is 0{^) and hence ^43 is negligible in the cases we shall consider. 

The reflected wave set up by an incident wave of unit amplitude and con- 
taining only the /•* mode (i.e. the incident wave described at the beginning 
of this section) is givenby the column matrix/" whose elements may be 
obtained from (3.3-12) and (3.3-13). Likewise, the transmitted wave is 
given by /+. 

PART IV 

GENTLE CIRCULAR BENDS IN RECTANGULAR WAVE 

GUIDES 

4.1 Propagation of Dominant Mode in a Gentle Bend — E in Plane of Betid 

When the magnetic intensity H lies in the plane of the bend, Hy = 0,' 
and equations (1.2-5) show that 5 = 0.. Thus we have to deal only with 
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A. In order to study the dominant mode we set / = in the cos {irty/b) 
{A depends on y through this factor) in the formulas of Section 1.3 which 
involve A and assume the dimensions of the guide to be such that h < a. 
We wish to determine -n, the first latent root of r„ defined by (1.3-5), 
from the approximate formula (3.2-2). In our case the elements 5m of diag- 
onal matrix To are obtained by putting «(= I) to zero in (1.1-5) : 

hi = rL = (7^ + {irm/a)\ ' m= 1,2,3, ■-• (4.1-1) 

so that (3.2-2) becomes 

yl := Tl + Fn - T,F,^F^,a\-\m' ~ iT' (4.1-2) 

The first task is to find the elements of the matrix F where, from (3.1-1) 
and (1.3-5), 

F=vl-Vl== (P-' - I)tI + P-'S (4.1-3) 

In the case under consideration P = / + J2 where J? is a square matrix 
whose elements are very small. In fact, the asymptotic expressions lead- 
ing to (Al-18) show that Ra and Sa are 0(^^), with ^ = a/pi , while i?,,- 
and Sij are 0(f) if J + i is odd and 0(?^) if i -|- _;' is even. When the approxi- 
mate value of P~^ obtained from (A2-2) is set in (4.1-3) and the matrix 
multiplications carried out it is found that 

/ °° \ °° (4 1-4) 

Fii = (-Rii+'L Ri^ R^i ) T% -hSu- E Rin. S^i + Oik') 

\ m=.l / m=l 

The "order of" symbol 0( ) will be omitted in the following equations, it 
being understood that the terms in the principal diagonal are correct to 
withm 0(|2) and the others to within 0(|). 

The values of the F's which enter (4.1-2) may be computed from the 
asymptotic expressions (Al-18) for the i^'s and 5's. They turn out to be 

Fir. = -4£m[4r?oT"'(m' - 1)"' + 3ar\rn - l)"'] 

. Pmi = -Hm[4Tl,^~\m' - ir' + a-\m' - l)"^] (4.1-5) 

Fn = e'[rL(l - 6ir-') + 60-V12 

In the expressions for Fi^ and F^i m is supposed to have the values 2, 4, 6, 
• • ■ . For odd values of m Fi^ and F^i are O(f^). When i = 1 in the ex- 
pression (4.1-4) for Fn, the two series therein reduce to Sz and Si where 

5, = E m{m' ~ IT" (4.1-6) 
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By expa'ndmg the typical terms in partial fractions and using the fact that 
the sums of (see, for example, page 238 of Reference^") 



(4.1-7) 



(4.1-8) 



for 5=2,4, 6, are 7rV8, 7rV96, 7rV960, it may be shown that 

S, = 7r764, S, - 7rV768 - 7rVl28, 

56 = (157r' - x')/3072. 

In (4.1-7) Bq denotes the q'^ Bernoulli number. The values of Sp may also 
be computed in succession from the two relations* 

Uip = Z-f 2 ' Cp4t-l,2i-l5p+f — Zj 2 Cp+t-l,2i-2 Jp4« 

where Cm.r^ is a binomial coefficient. Still another method is to make 
use of the generating function 

E ^5p+i = (1 + /) S (*»' - 1 ~ ^)~' = 1 - hT^^ cot 7r.i; 

p=,0 jii=2,4,e--- 

where 4a;^ = 1 + ^ Note that by this definition Si is J in contrast to the 
non-convergent series obtained by putting ^ = 1 in (4.1-6). 

Substituting the values for the F's given by (4.1-5) in the expression 
(4.1-2) for 71 and using the sums (4.1-8) of the series which occur gives 

J[ = rl, - -4 [1 + aV'M - 6t-') + (ario/,r)\5 - t73)1 (4.1-9) 

When the dominant mode is propagated without attenuation both 71 and 
Tio are negative. 

The general form of (4.1-9) has been obtained by both Euchholz^ and 
Marshak^ by different methods. In our notation their result is 

tL = rL - 5 [1 + .'rid - 6.-'»r') + (^")' (5 - .'»V3)] 

(4.1-10) 

where 7„i„ is the propagation constant for the m, w*"" mode when the mag- 
netic vector is in the plane of the bend. 

* I am indebted to John Giordan for these relations. 
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4.2 ReHection Due io Dominant Mode Incident upon Gentle Bend — H in 
Plane of Bend 

Let the system be the one described in the first paragraph of Section 2.3 
and let the incident wave contain only the dominant mode. Then the 
matrix propagation constant is the T^ of Section 4.1 and the column matrix 
h specifying the incident wave has unity for its top element and zero for its 
remaining elements, i.e., /> = 1 in the formulas of Section 3.3. 

We shall be interested only in the reflection coefficient, jT of the dominant 
mode. Here we shall denote it by gi^,, in line with the notation of equa- 
tion (1.1-3), in order to distinguish it from the corresponding coefficient 
(which will be denoted by rfoi) when E lies in the plane of the bend. ,jV 

Setting ^ = 1 in the expression (3.3-13) for the reflection coefficient and 
using equation (4.1-1) for 5™ gives 

■ /T = gro = -e'^*''""'''Ui(sinh2c7i)/2 +..42] ' (4.2-1) 

where 71 has just been obtained in (4.1-9) and . , , 

Ai = 2iiii + (7? - r?o)rro' 



m=i L T^a \m^ — 1) J 



^. = E 



cosh 2c7i - e ^'■*''' 



(4.2-2) 



ti2VmVioTr^a-\m? -1) 

• {VlmFml^mQ + VmlFlm^Vi + FlmFnii\ 

From (Al-18) and Fu = 1 + vn it follows that 

Vrm= 7^1 = %T-'kmW ~ 1)-=' 



(4.2-3) 



where m = 2, 4, 6, ■ ■ • . For odd values of m, Vi^ and F„,i are 0($^)* 
Substituting these values together with those for the F's given by (4.1-5)' 
using the sums (4.1-8) and the expression (4.1-9") for 71 — Fio finally leads 
to (after considerable cancellation) 

A,= - ?'rroV74 ; '■' (4.2-4) 

Likewise, for even values of m, 

VimFrnxVL + V^iFiA + FimFmi = leem'a-^m^ ~ 1)-^ (4.2-5) 

All of the terms in the expression (4.2-1) for ^ro are now known (the values 
of 7m may be obtained by setting m = in (4.1-10)). We shall make the 
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further approximation of putting r„,o for 7„ . Since r„o — y„, is 0($2) no 
serious error is introduced and we have 

_ _ g^ sinh 2crio 1^ ^ (cosh2crio-e-='''"°) tn . . 

gio — (jp2" 2 — a Zj r V ^"^ ri TTa ^4.2-6} 

Ol loa IT m=2.4.S-.. ilO 1 ...0 O (w^ _ 1J3 

in which 

aT^o = [T-(m^ - 1) + ahio]', ^ = a/pi . (4.2-7) 

- For frequencies such that only the dominant mode is propagated the 
ratio of the power in the reflected wave to the power in the incident wave is 
I gTo p. Marshak has given an expression for this ratio which is the same 
as that obtained from (4,2-6) when the neghgible (for his case) terms e"^'^ ^'"° 
are omitted. 

The corresponding expression for the transmission coefficient derived 
from (3.3-14) ior ft is not as simple as (4.2-6). 

4.3 Propagation of Dominant Mode in a Gentle Bend — E in Plane of Bend 

When the electric intensity E lies in the plane of the bend, Ey = 0, and 
equations (1.2-5) show, that A = 0. Here we deal with B in much the 
same way as we dealt with A in Section 4.1. The dominant mode is ob- 
tained by setting / = 1 in the sin{irfy/b) in the formulas pertaining to B 
in Section 1.3, It is assumed that b > a. 

Examination of the matrices (1.3-14) indicates that, for the sake of con- 
venience, we should call the top row of our matrices the O' row and the left- 
most column the column. In line with this we call 70 the propagation 
constant of the dominant mode in the bend. The elements 5,^, of the diagonal 
matrix To are obtained by putting n{= /) = 1 in (1.1-5): 

bl = rL = tr' -I- {mn/af -\- -K^lh"-, w = 0, 1, 2, • ■ • (4.3-1) 

When we make the appropriate shift in the subscripts, equation (3.2-2) 
yields 

CO 

70 = To, -I- /'oo + E /^o„. /^,„o a' tt"' m"' (4.3-2) 

in which the elements of the matrix ¥ are to be determined from (1.3-13); 
F =v\-~ V\ = {Q-' - I)tI + Q-'U (4.3-3) 

As in (4.1^) we have, with Q = I + T, 
Fii= -TiiV),-{- Uu 

(4.3-4) 



\ m=0 / m=0 



u„ 
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By using the asymptotic expressions (Al-19) for the T's and U*s and sum- 
ming the series with the value of (4.1-7) for g =4 given in Section 4.1 we 
obtain 

F,,o = -S^tIt-'m-^ (4.3-5) 

In these exjM-essions m is supposed to have the values 1, 3, 5, • ■ ■ . When 
m is even Fom and /^,„o are 0(?^). 

Substituting (4.3-5) in (4.3-2) and summing the series with the heW 
of the values of (4.1-7) given in Section 4.1 gives 

70 = tI - ^'tI{5 + 2a'vl)/60 ' (4.3-61 

A result equivalent to (4.3-6) has been given by Buchholz who also gives 
the approximation to the propagation constant when m >0 (and the elec- 
tric vector in the plane of the bend). In our notation his approximation 
is 



7™n im« ~^4a''L 



\ 7rm / 

'or, 



+ i 



TTtn 



y(21+7r^m^)] 



(4.3-71 



In writing (4.3-7) we have corrected a misprint in Buchholz's expression. 
In order to agree with Buchholz's equation (5.30a) the leading term within 
the square brackets would have to be changed from 3 to —3. This change 
was mdicated by the results obtained when our equation (3.2-2) was used 
to obtain special cases of (4.3-7). Probably the best way of obtaining 
(4.3-7) is furnished by Marshak's method (WKB approximation, out to 
second order terms, applied to BesseFs differential equation). If one wishes 
to verify (4.3-7) by using Marshak's report^ as a guide, he should correct 
the misprint in Marshak's equation (12a). 

4.4 Reflection Due to Dominant Mode Incident upon Gentle Bend — E in Plane 
of Bend 

The problem here is the same as that treated m Section 4.2 except that 
now the electric vector hes in the plane of the bend. In line with equation 
(1.1-4), the reflection coefficient /r for the dominant mode will be denoted 
by d^i . As in Section 4.3 the subscripts indicating the position of matrk 
elements will be adjusted so as to start with instead of 1. The square 
matrix W given by (1.4-5), and associated with the junction conditions for 
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B in the same manner as V is associated with A, now replaces V. Thus 
our expression (3.3-13) for the reflection coefficient becomes 

/o"i = (^ = - [^i(sinh 2cT,,)/2 + A2] (4.4-1) 

where we have neglected the difference in Foi and 70 and where 

Ai = 2wm + (70 - Tli)Toi - Z ■ 

n om n',„o -r — T-^ — ^ 



(4.4-2) 



^2 = z 



cosh 2croi — e 



-2c r. 



1 2r„iroix''a-2m2 

From 1^00 = 1 + w^oo and the asymptotic expressions (Al-19) it follows 
that, for m = 1, 3, 5, • • • 

Woo = ^Vl2 

(4.4-3) 

For even values of m, W^m and H-^^^o are 0(^^). Substitution of these values 
together with those for the T^'s given by (4.3-5), using the sums (4.1-7) 
and expression (4.3-6) for 70 — Fqi leads to 

^1 = ^V12 

Tro».F„,orL + W„^,F^vli + F^„J,„^ = - U\l,ar\~''m-^ 

for m odd. 

Thus the reflection coefficient for the dominant mode when E lies in the 
plane of a gentle bend of length 2c is approximately 

^^sinh 2croi -, ^^4roi v" cosh 2croi - e"^-'^"^ ,, , ,, 
d^i = -^ + — -^ 2^ -T^ (4.4-4) 

where Tmi is given by (4.3-1) and b > a. 

PART V 

NUMERICAL CALCULATIONS 

5.1 Bend in Plane of Magnetic Vector 

Let a/b = 2.25 and Xo/a = 1.400 where Xo is the free-space wavelength 
of the dominant wave striking the bend. The propagation constant 
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Tio of the dominant wave is obtained by setting m = 1, » = in (1.1-5). 
The r's corresponding to the higher modes may be obtained from the 
same formula: 

fl^rL = - (^Y + TT- = -10.272 aV,, = / 3.205 

aVL = a'r'io + 3/ = 19.336 aV^^ = 4.397 (5.1-1) 

0Y30 = o'r?o + 8/ - 68.684 aFso = 8.288 

We shall consider a 90° bend. The approximation (4.2-6) appropriate 
to gentle bends becomes 

g- = if [- .0122 sin (5.03/?) + .0087 cos (5.03/?)] (5.1-2) 

where the exponential terms have been omitted since they are generally 
negligible. In (5.1-2), ? = a/pi and the arguments of the sine and cosine 
terms arise from IcVia = ■jrario/(2?). From (4.1-9) the approximate 
change in the propagation constant produced by the curvature is obtainable 
from 

7? - r?o = .294?Vfl' (5.1-3) 

where 71 is the propagation constant of the dominant mode in the bend. 

The determination of gfo by matrix methods will be illustrated for a 90" 
bend in which pi/o = 0.6. This makes c/a = pi5r/(4o) = .4712, crio = 
tl.510 and the appropriate equations in (2.3-5'l and (2.3-4) become, upon 

setting /r = i?ro and /i = gio , 

g+ = e-'^'^xi - yi) = (.061 -{- j.998)(:ci - yO 

g- = /^'»(j:i + yO - /='- = (.061 + i.998)(:ri + y,) (5.1^) 

+ .993 - J.121 

Here Xi , yi are the top elements in the column matrices x, y. The problem 
is to compute x and y from the matrix equations (2.3-3) with T replaced 
by Ta , To defined by (1.3-2) with ^ = 0, and k a column matrix whose 
elements are zero except the top one which is unity. Since the order of the 
matrices is infinite, an exact solution calls for an infinite amount of work. 
A compromise must be made between the accuracy desired and amount of 
work one is willing to do. The following numerical work uses third order 
matrices. 

The first step is to compute the square matrix, obtained from (1.3-5), 

a'rl = P~\a-Tl + aS) (5.1-5) 
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The elements of the diagonal matrix a^rl are given by (5.1-1) and those of 
P and S by the equations and tables of Appendix I. 

3.3447 1.6087" 



aTi = 



1.4429 .8812 .6821 

.8812 2.1250 1.3745 

. .6821 1.3745 2.4879 

■-9.086 -1.785 

.157 17.218 

.996 -13.566 



-12.292 
-6.3039 16.369 6.9738 
-3.5034 -11.3031 65.213 
-5.178" 
-19.362 
38.329 



-^ (S.1-6) 



The next step is to use (5.1-6) to evaluate the coefiicients of x and y in 
(2.3-3). The square matrices T^c tanh T^c and r„c coth V^c cause most of 
the computational difficulties. We shall deal with these matrices by using 
Sylvester's theorem (an account of this theorem is given in Section 3.9 of 
Reference®). This requires the determination of the latent roots and modal 
rows of aV],. However, it is interesting to note that the matrices in ques- 
tion may also be computed from cT\ (which is easily obtained from aVt) 
by processes which employ only matrix multiplication, addition, and in- 
version. 
Thus, setting A^ for c^Ta, 

^->sinh>l = ^ +^ + 57+ •• 

cosh ■^=^+yj'"^4j + '' 
A coth A = (cosh A)(A-^ sinh /()- ' 
A tanh A = A^{A coth A)-\ 

Although the series always converge, they do so too slowly to be of use in 
our computations. The same is true of the series 

A tanh X = E 8^- [{2m - D'ttV + 4^^^]-'. 

For the matrices we shall encounter it appears best to use Sylvester's 
theorem even though this requires the determination of the latent roots and 
modal rows of aT\ . The square matrix formed from the modal rows* 
will be denoted by n. 

* As has already been mentioned in the footnote associated with equation (1.3-9), we 
shall use the notation and theory set forth in Sections 3.5 and 3.6 of Reference*. 
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We shall use the relations* 

VaC tanh V^c ^ «"^' [7c tanh 7c]d k 

r„c cosh PqC ^ K^^ [7c coth 7c]d k 

where the subscript d on the brackets stands for "diagonal" matrix, the 
i**" element in the principal diagonal of [7c tanh 7£:]d being 7 iC tanh7iC where 

7.- c = c\T/a (5.1-8) 

and Xi is the t'th latent root of F^ . In our applications 7^ is either posi- 
tive real or positive imaginary. 

From (5.1-6) the Xj's are the roots of 

85 5.178 

(5.1-9) 



X + 9.086 




1.785 




5.178 


-.157 


X - 


- 17.218 




19.362 


^.996 




13.566 


X - 


- 38.329 



- X^ - 46.461 X^ - 101 .96X + 3464.5 = 

and have the values 

Xi = -8.886, X2 = 8.284, X3 = 47.06 (5.1-10) 

The elements K21 , Kai of the modal row [1, khi , K31] corresponding to X: 
may be obtained by solving the two equations derived from the last two 
elements of 

[\, K2, , K^,]{\J - a'vl) = (5.1-11) 

namely, 

1.785 + (Xi - 17.218) K21 + 13.566 K31 = 

5.178 + 19.362 K21 + (X: - 38.329) K31 = 

When the value of Xi from (5.1-10) is used these equations yield 

K21 = .1593, K31 = .1750 

Likewise, the first and third elements of 

[ku , 1, 'C32](X2/ - aVl) = 

and the first and second elements of 

[ku , «23 , IKXa/ - arl) =0 

*This is the modal row matrix analogue of equation (11) in Section 3.6 of the Refer- 
ence^. The modal rows of Tn are equal to the modal rows of a^r^o. 
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give , - 

KI2 = .0465 K32 = .6524 

Ki3 = .0165 K23 = -.4555 

Thus, the numbers entering (5.1-7) are 

7ic= .4712 (-8.886)'" -n.404, -y.c = .4712 (8.284)"' = 1.356 



7ic tanh 7ic = —8.382 
Yiccolh'^ic = .2354 



ync tanh 72C = 1.187 
72C coth 72C = 1.549 



7ac = .4712 (47.06)"' -- 


= 3.233 




"1 


.1593 


.1750 


73C tanh73C = 3.228 




K = 


.0465 


1 


.6524 


73C coth 73C = 3.243 






_ .0165 


— .4555 


1 



For the purpose of calculation it is convenient to transform (2.3-3) by 
inserting (5.1-7) and premultiplying by kF^\ We obtain 

([7c tanh yc\dii + k\^^\\c) x = KV~'^cV^i "h 

([7c coth yc]dK + kF"Voc) y = Kl^'croe'"'"// 



(5.1-12) 



in which 



kV-' = 



1 

.0465 
.0165 


.1593 
1 
- .4555 


.1750 
.6524 
1 




1.1204 .3911 .1629 

.3911 1.2833 .4946 

_ .1629 .4946 1.3460 


.9492 
- .2608 

.1427 


-.1992 

.7686 
- .7900 


.0883" 

.2339 

1.0160_ 


! 


Tioc - i 1.510 
r2oC - 2.072 
Taoc - 3.905 



where the elements of V are obtained from the formulas and tables of Ap- 
pendix I. 

The i^ equation of the set obtained by writing out the first of equations 

(5.1-12) is 

E Uii-yiC tanh y^c + {kV-^) i ,V ^.c\x ^ = {KV-')ncVa-''''' (5.1-13) 
'■=1 ! ; ' ■}.'. j' 

where {kV^^)^ denotes the element in the i*'' row and/^ column of kV~\ 
K,i is the element in the i^'^ row and /^ column (note the reversal of the 
usual convention regarding the order of subscripts) of k, kh = 1, and h 
has disappeared because it is a column matrix whose top element is unity 
while the remaining elements are zero. It will be noted that the only 
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imaginary terms in (5.1-13) occur i;^ the coefficients of xa and arise from tlie 
imaginary quantity Tioc. By making the substitution 

'^' l+«i.rio ^^-^ ^^ 

the set (5.1-13) may be reduced to 

((ciiTiCtanh yic)ui + 2^ [KayiC tanh^iC + {K.V^^)ijVj^c\ui = (kF^^),-! 

(5.1-15) 

in which the coefficients are all real. It should be noticed, however, that 
nothing is gained by making the substitution (5.1-14) when the frequency 
is so high that other modes in addition to the dominant are propagated. 

The equation for y corresponding to (5.1-15) may be obtained by re- 
placing tanh by coth and uhy v where now 

Incidentally, if we set j = 1 in (5.1-14) and (5.1-16) and substitute in 
the expressions (5.1-4) for gio we may show that, since u-y and vi are real, 

Uu.r + krof = 1 (5.1-17) 

Equation (5.1-17) may be obtained at once from the fact that the energy 
of the waves leaving the bend must equal the energy of the incident wave. 
It may also be shown that ^ro vanishes when uiVic^Via = 1. 

When the above numbers are set in the three equations obtained from 
(5.1-15) we get 

-8.382 Ml -1.748 mj -1.122 u^ = .9492 
.055 Ml +2.780 «2 -1-1.688 M3 = -.2608 
.053 Ml -3.105 M2 +7.191 uz = .1427 

from which 

uy= - .0940, xi = .1400 + i.0113 

The equations for ^1 obtained by substituting coth for tanh are 

.2354 vi - .3732 v^ + .3861 V3 = .9492 
.0717 vi + 3.1417 V2 + 1.924 vz = - .2608 
.0534 vi - 3.1143 V2 + 7.211 ^3 = .1427 



from which 



vi = 3.930, yi = -.1045 + j.9803 
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When these values are set in (5.1-4) we finally obtain 

gto = .9822 + T.1858 
gw = .0048 ~ 1,0255 

The following table lists values ofgfo obtained by the methods of this 
section*. Here the bend is in the plane of H, a/b = 2.25, Xo/a = 1.400, 
where Xo is the free space wavelength of the incident wave, pi is the 
radius of curvature of the axis of the guide. The smallest possible value of 
pi/a is 0.5. The term "approx." refers to equations (5.1-2) while "1st 
order", "2nd order", etc. refers to the order of the matrices used in the 
computations. The amplitude of the reflected wave is gTo and the ampli- 
tude of the wave sent forward is gto when the incident wave is of unit ampli- 
tude. 



Pi/a Approx. 

.6 

.7 

.8 

.9 
1.0 
1.2 
1.5 

.6 -1.0280 .0020 -i.0074 .0056 -1.0280 .0048 -t.0255 

.7 -i.0068 -.0005 +J.0023 .0013 -f0131 .0007 -1.0066 

.8 +1.0062 -.0013 +i.0074 -.0003 +i.0039 -.0004+1.0051 

.9 +2.0128 -.0014+1.0087 -.0009 +i.0123 -.0009+1.0123 

1.0 +i.0143 -.OOlO+i.0075 -.0010 +i.0148 -.0010 +i.0147 

1.2 +1.0079 -.0002 +1.0018 -.0005 +t'.0O86 -.0005 +i.0085 

1.5 -1.0040 +.0003 -1.0034 +.0002 -1.0041 +.0002 -1.0042 

It appears that the values obtained from the first order matrices are quite 
far from the true values. On the other hand there is considerable agreement 
between the approximation and the second and third order values, especially 

at the larger values of pi/a. 

5.2 Bend in Plane of Electric Vector 

The calculations for this case are quite similar to those presented in Sec- 
tion 5.1. If we are to deal with the same waveguide it is necessary to 

* The computations were performed by Miss M. Darville. I am also indebted to her 
for the values given in the tables in Section 5.2 and Appendix I. 





gio 










1st order 


2nd order 


3rd order 


.964 


+i.267 


.980 


+1.197 


.982 


+i.l86 


.974 


+i'.224 


.994 


+1.105 


.994 


+1.111 


.984 


+1.178 


.997 


+1.082 


.997 


+1.082 


.988 


+1.153 


.997 


+1.074 


.997 


+1.073 


.991 


+1.135 


.998 


+1.066 


.998 


+1.066 


.994 


+1.110 


.998 


+i.056 


.998 


+1.056 


.996 


+t.084 


.999 


+1.043 


.999 


+1.044 
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interchange the dimensions a and h so that now h/a = 2.25 and, for the 
same frequency, \<i/b = 1.400. 

For a 90° bend the approximation (4.4-4) for the reflection coefficient 
d^i for gentle bends, i.e. for ? = o/pi small, becomes 

' d~,i = ?^2[-.04l7 sin (2.23/^) + .0209 cos (2.23/^)] 

where the negligible exponential terms have been neglected just as in the 
analogue (5.1-2) for gTo ■ 

The following table, which is sunilar to the table at the end of Section 
S.l, gives the results of computations for bends in the plane of the electric 
vector. 

Pi/a Approx. 1st order 2nd order 

.6 .823 -\-i.547 .975 -\-i.223 

.7 .887 4-i.447 .994 +i.051 

.8 .921 +J.380 .996 +i.042 

.9 .941 +i.332 .997 +1.035 

1.0 .954 +i.295 .998 +1-031 

1.2 .970 +i.242 .999 +t.023 

1.5 .982 +i.l90 1.000 +i.017 

.6 -i.0996 - .0855 +t.l284 - .0020 +i.0086 

.7 -i. 0848 - . 0520 +i . 1031 + . 0050 - i . 0975 

.8 -i.0706 -.0330 +j. 0800 +.0033 -i.0792 

.9 -t.0575 -.0214 +i. 0605 + .0022 -t. 0635 

1.0 -i.0457 -.0137 +i. 0443 .0021 -i.0507 

1.2 -t.0258 -.0051 +i. 0204 .0007-1.0282 

1.5 -i.0051 +.0001 -1.0004 .0001 -i.0062 

The agreement between the approxunation for d^ and its second order 
matrix value is fairly good from pi/a = .7 onward. 

APPENDIX I 
Calculation of Pp,,, etc. foe Circular Bend 
It is convenient to write Pp,n and Qp,., as given by (1.2-10) and (1.2-15) 
in the form -■ '-'^ ■ 

Pp^ = 5^ + -Rp™., P,ffi = 1,2,3, •v:- .. ■ 

- , ■ ■ (Al-lJ 

Qk«. = hl+T^, p,m = Q,i,2, ■•• - ■ • 

where Sm is unity if ^ = w and is zero otherwise and .1 ,.; ■ . 



(Al-4} 



REFLECTIONS FROM CIRCULAR BENDS 341 

Rpm = (2/a) I (pip"' - 1) sin (x/'.v/o) sin (jrm.r/o) dx 

(Al-2) 

Tp,„ = (€p/a) j (plp~' - 1) cos (irpx/a) cos (irmx/a) dx 
Jo 

in which cd = 1, e? ^ 2, ^ = 1, 2, ■ ■ ■ 

In (Al-2), (1.2-11), (1.2-16), (1.4-3), (1.4-6) we make the substitutions 

P — m = r, u = pi tt/o — 7r/2 ^ pax/o 

P -\- VI = s, V = pi ir/a + t/2 ^ Pst/*^ (Al-3) 

y = TTx/a, w = piT^/o, p = .V + pi — a/2 = a{y + j()/7r 

Introduction of the integrals 

h r (lA) / [w'(>' + »)"- - 1] cos .?>< dy 
Jo 

= (1/a) / (pip"- — 1) cos (TTSx/a) dx 

J, = T I -^-^ dy = TT sin (wsx/a) dx/p, 
Jq y -\- ii Jo 

a Jo y -r It 

enables us to write 

R„. = Ir- L, S,,., - -ma-%r, + /,) 

T„„, = £,(/. + Q/2, Uy,n = me^a-KJ, - Jr)/2 (Al-5) 

Kp™ = Kr - K, , Wp„^ = ep{Kr + K,)/2 

where /^ and K, are even functions of .? and /, is an odd function of j. 
£0=1 and Cp = 2, ^ = 1, 2, 3, ■ ■ ■ . Since w and « depend only upon the 
ratio pi/a, the values of I, , K, and /, depend only upon pi/o and the 
integer s. These quantities are tabulated at the end of this appendbc. 
Setting y -\- u equal to t gives 

/, = TT f sin s(t. - u) dt/t 
J It 

= Tr[Si{sv) — Siisu)] cos su - Tr[Ci(sv) — Ci{sii}] sin su 

where Si and Ci denote the integral sine and cosine functions. Integrating 
by parts enables us to express Is in terms of A . Thus 



Jq 



(y + ii)~' cos sy dy = tr'^ — v ^ cos sir — it ^ sJ^ (Al-7) 
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and 

I, = TT-'w'^br' - (-)Vi - T-'sJs] CAl-8) 

except when ^ = in which case 

/(, = TT-^wKtr' - v~^) - 1 = W^/(21V) - 1 

(Al-9) 

= 7rV(4M^) = [(2pi/ay - 1]-' 

When pi/a is large, u and v are large, and the asymptotic expansion of 

(A 1-6) gives 

w-'j. - s-'[u-' - (-)V'] - 2\s-\u-' - (-)V^l + -■■ (Al-10) 
When (Al-10) is placed in (Al-8) 

/, - ^-'2\wh-'[u~' - i-Yv-'] - ■■■ (Al-U) 

Formulas for Ks may be obtained in much the same way. 

K> = (pi/a) / cos5{/ -u) dt/l 

= {pi/a){[Ci{sv) - Ci(su)] cos su + [Si{sv) - 5i(sM)lsin su] 

and when s = 

K, = (pi/fl) log (1 + 7r/«) (Al-13) 

The asymptotic expression is 

ap-^'K. - s-\u-' - (-)V^] - 3!r*k"* - (-)"rl + • ■ ■ 
It is convenient to write the asymptotic expressions in terms of the new 
variable 

i - fl/pi (Al-14) 

When s is even and greater than zero 

/« - e/s, L -- 6fV-V^ K, ^ lei^-H-^ (Al-15) 
and when j is odd 

/, -- l^s, L — 4^x-V^ K, r^ 2^x-2j-2 (Al-16) 

When 5 = 

/o - ?V4, ffo - 1 + £V12 (Al-17) 

We shall need the following asymptotic expressions which may be obtained 
from the above work 
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5i„, S„,i , Si,, ■^4a-'^m{ni^ - Ij"' 

where m = 2, 4, 6, • • • . Also, if now m = 1, 3, 5, ■ ■ ■ , 

Too - ^V4, t/oo = 0, T'l^oD = 1 + ?V12 
r„o '^ 8^m-V-^ (7.0 = 
To™ ~ 4$m-V-^ Uon. ~ 2a-2^ 
I^on, ~ 2$m-V-^ lK„o ~ 4^m-^Tr~^ 
Values of la , Ja and Ka 



(Al-18) 



(Al-19) 





PiA 




h h 


h 


h 


/* 


h 


h 






,5 




so 00 


00 


CO 


00 


CO 


CO 






.6 


2 


.2723 2.31879 


1.82979 


1.43755 


1.14772 


.94432 . 


,78458 






.7 


1.04166 1.28232 


.76232 


.53315 


.37692 


.28832 . 


22090 






.8 




.64103 .88256 


.44628 


.29206 


. 18995 


.14101 . 


09986 






.9 




.44643 .68200 


.30111 


. 18992 


.11546 


.08517 , 


,05908 






.0 




.33333 .56052 


. 22008 


. 13637 


.07818 


.05814 . 


,03865 






.1 




.26041 .47844 


. 16916 


. 10459 


.05772 


.04298 , 


,02759 






.2 




.21008 .41905 


. 13486 


.08413 


.04344 


.03361 . 


,02056 






.3 




.17361 .37361 


.11059 


.06961 


.03500 


.02732 . 


,01633 






.4 




.14620 .33780 


.09259 


.05926 


.03168 


.02293 , 


,01297 






.5 




.12500 .30876 


.07872 


.05147 


.02315 


.01969 . 


,01068 




2.0 




.06667 .21803 


.04133 


.03080 


.01135 








2.5 




.04167 .16980 


.02564 


.02217 


.00675 












/o 


/i 


/. 


/i 


/* 


/i 


/t 




.5 





















.6 







3.99809 2.01979 2 


.31576 1 


.48355 


1.66348 


1.15716 




.7 







3.21624 1. 


3054 1 


.58176 


.84936 


1.04899 


.61931 




.8 







2.72356 . 


73339 1 


.21541 


.56683 


.77644 


.40134 




.9 







2.37231 . 


70698 


.99327 


.41079 


.62183 


.28546 




.0 







2.10615 


55663 


.84343 


.31379 


.52170 


.21578 




.1 







1.89624 . 


45091 


.73508 


.24849 


.45123 


. 16981 




.2 







1.72581 . 


37335 


.65279 


.20254 


.39869 


. 13768 




.3 







1.58448 . 


31450 


.58812 


.16843 


.35788 


.11413 




.4 







1.46507' . 


26878 


.53573 


.14216 


.32515 


.09636 




.5 







1.3628 


23251 


.49238 


.12243 


.29825 


.08254 


2 


.0 







1.0122 


12817 


.35299 


.06596 






2 


.5 







.8062 


08128 


.27660 


.04140 
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Pi/a Ko Ki Ki Ki Kx K^ Ke 

',('. .5 00 00 00 00 OO 00 CO 

.6 1.43874 .61659 .31832 .22547 .15539 .12199 .09275 

.7 1.25423 .42303 .17355 .11508 .06971 .05353 .03704 

.8 1.17306 .33141 .11439 .07463 .04093 .03195 .02038 

.9 1.12748 .27575 .08261 .05431 .02737 .02214 .01325 

1.0 1.09861 .23761 .06305 .04244 .01976 .01675 .00938 

1.1 1.07891 .20953 .04994 .03471 .01509 .01341 .00705 

1.2 1.06476 .18785 .04075 .02935 .01192 .01126 .00548 

1.3 1.05421 .17050 .03391 .02539 .00977 .00956 .00443 

1.4 1.04610 .15626 .02871 .02243 .00807 .00838 .00365 

1.5 1.03972 .14434 .02467 .02013 .00682 .00745 .00315 
2.0 1.02165 .10508 .01332 .01330 .00358 

2.5 1.01366 .08295 .00838 .01010 .00217 

APPENDIX II 

Functions of Almost Diagonal Matrices 

Let £ be a matrix whose elements are small in comparison with unity. 
It is then often possible to approximate a matrix deiined as some function of 
the matrix I -\- E, where / is the unit matrix, by the expansion 

/(/ + £) = 1/(1) + ^/'(l) + ^V"(l) + ■ • ■ ■ (A2-1) 

Thus, for example, when we take/(s') to be z~^ we obtain 

(/ + E)-^ ^ I ~ E + E"- ■■■ . CA2-2) 

Here we shall give similar formal results ior f{D + E) where now D is a 
diagonal matrix 



D = 



di ■ 
d2 ■ 

dn 



(A2-3) 



whose diagonal elements are unequal and the elements £<, and Eji are small 
in comparison with the absolute value of 1 t^i — (/, | . We shall restrict 
ourselves to a first approximation of the non-diagonal terms oi f{D -\- E) 
and to a second approximation of the diagonal terms. The results are closely 
related to the ones obtained from the perturbation theory used in wave 
mechanics. 
We assume that/(Z) + E) may be defined by the series 

/(D -VE) = a,I + oi(Z) + £) + a2{D + £)>=+••• (A2-4) 
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where o„ is a scalar and 

(D + Ey ^ (D-{- E)(D -^ E) = D^ -\- DB -h ED + B? 

and so on. The sum of the terms independent of E is /(£>). The terms 
of order E are 

EmD-\~ E 

DE + ED in (D + Ef 

(A2-5) 
Dm + DED + ED""' in {D •]- Ef 

2 D^ED"" in (Z> + £)" 

where the summation extends over the non-negative integer values of t 
and m for which l-^- m = n— 1. The element in the ith row andj'th column 

f f 

of D EJT is diEijdJ and hence the corresponding element in the summa- 
tion in (A2-5) is 



(di — d")/(di — dj), ij-ij 



EiM^ = „_i (A2-6) 

[ndi , i = J. 

Thus the terms of order E in the ith row and jth column of f{D + E) are, 

from (A2-6) and (A2-4), 

'ii - <i7 (A2-7) 

Euf (di), i ^ j 

where the prime on/ denotes its first derivative. 
The terms of order E^ in (Z) + £)" are 

E d''ED^ED'^ - E [dtEiMd^iEiid^U (A2-8) 

where the summations extend over all the non-negative integer values of 
k, I, m for which k \ t -\- m = n — 1. On the right [(fi£,-,-]M denotes a 
square matrix whose element in the ith row and ^th column is diEij . 
Likewise the second factor in brackets is a matrix having diE^S] in the ith 
row andj'th column. The element in the ith row and jth column of (A2-8) 
is, from the rule for the product of two matrices, 

i; E (rf'£i,)((/f£.y^) = E £.>£.; E d\d[ri. 



If i, 5, andj are unequal the sum in U, t, m is 

1 pr - d: _ ^ ■ - dn 

di — dj [_ di — da dj — ds J 
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and in case of equality the sum may be found by a limiting process. Since 
we are interested in this order of approximation only for the diagonal terms 
we set i = j and obtain for the sum 

di - d, (di - d,y' 
n(n ~ 1) „_3 

Thus the contribution to the i\h diagonal element oif{D + E) from terms 
of type (A2-8) is . , 

where the prime on S indicates that the term s = ? is to be omitted. 

Thus, to summarize, we may say that the first approximation to the non- 
diagonal term in the ith. row and^th column (i 9^ j) oif{D + E) is 

di — fly 

and the second approximation to the diagonal term in the ith row and Jth 
column of /(Z) -\- E) is . ■ . . 1 

m+E,,f'(d) + ^f"{d,) 

(A2-11) 






where the primes on / denote derivatives and the prime on Z indicates that 
the term ^ = j is to be omitted. 

Two results obtained from (A2-10) and (A2-11) are of interest. For 
the first result we set/{2) = z~^ and get the following approximations to the 
elements of (D + E)~^: 

~ Ei,{dai)-\ i^j 

2^ ^ 1I <A2-I2) 

d7 — d7 \ Eii - 2-. EuE,,d^ , i = j. 

For the second result we set /(s) = s ^ and obtain the following approxi- 
mations to the elements of {D + £)"^: 

E,MT + df)-\ i^j - 

dr + K"^ [s. - 1 E^.E^^d'r + dTr'\ , i = j. ^^^"^^^ 

In (A2-12) and (A2-13) the summations include the term s = i. 
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We shall now state several results related to the above formulas. Let w 
denote the matrix Z) + £) so that the tj'pical element w,/ = £,-; , i ^ j, 
and ;*ii = di-\- En . Then the latent roots Xi , As , - - ■ Xjv of m are the roots 
of the equation obtained by setting the determinant of X/ — « to zero: 

[X — »11 —Ul2 

I X/ — « I = i — U-:i X — "33 



= 0. 



(A2-14) 



The modal column k,- corresponding to the jth root X,- satisfies the matrix 
equation 

iXjT - u)ki = 0. (A2-15) 

Since the non-diagonal elements of ii are small, we see from (A2-14) that 
we may label the roots so as to make X; nearly equal to itjj , and this together 
with (A2-15) shows that all the elements of fe,are nearly zero except the_;th 
which we may choose to be unity. When these approximate values are 
taken as a first approximation in the process of solving {A2-15) by successive 
approximations, the second approximation is found to be 

hi 



ki = 



_RNj_ 



(A2-16) 



k ■ = 



Ui, — Us 



s^j 



where the 1 in the column for kj occurs as the jth element. This expression 
for Xy occurs in the perturbation method often used in wave mechanics. 
For the modal row k,- corresponding to X; we have in much the same way 

Ky(Xy/ - u) - 

Ki ^ [Kij , Kij , • ■ ■ Kjj , • ■ ■ KNj] (A2-17) 

Ujs Kjj 



Ksi 



Uji — 1t,s 

where the last expression is an approximation and where kj, may be chosen 
at our convenience. 

The results (A2-10) and {A2-11) may also be obtained from (A2-2), 
(A2-16) and the relation* 

7{Xi) ■ 

/-(Xa) - 



/{«) = k 











/(x^). 



(A2-18) 



* This is equation (12) in SecLion 3.6 of Referenced Although proved only for poly- 
nomials it may be verified to be true for the applications which we shall make. 
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where k is the square matrix [ki ,ks , ■ ■ ■ kn] composed of the modal columns 
kj , and the latent roots are assumed to be distinct. The element in the ith 
row andy column of k~^ is, to our order of approximation, equal to —kij 
when i 7^ j and to 

1 +E'^,a«; (A2-19) 

8=1 

when i = j. Here, as usual, the prime on S indicates that the term s = j 
is to be omitted. 

Another result which is sometimes useful is the approxunate solution of 

UiiXi + Mia-Va + ■ ■ ■ + UipXp + ■ ■ ■ niKX?/ = tji 

Jfai^'i + U22X2 + • • ■ + M2j,a;p + ■ ■ ■ + U2nXk = V2 (A2-20) 



UpiXi + Up2X2 + ■ ■ • + UppXp + ■ ■ ■ + tlpffXfi = ?;, 



where \vi\^\vp\J ^ P) ^.nd the m's are of the type assumed in (A2-14) : 

x,=^ - !!il^^ j ^ p (A2-21) 

If 
a:p = — + S' — ^ (Vpi*,p - i?«Mpp). 

upp »=1 Wgg Mpp 

The prime on S indicates that the term s = p (which happens to be zero 
here) is to be omitted. A somewhat more general result which states that 
the ith element in the column matrix (D + E)~ t] is 

Vi d7' +J1 Ei. dT' (r;\E.i vi - di 7,,) {A2-22) 

may be derived from (A2-12). Here the elements of the column matrix 
need not be restricted by \-ni\ '^\'np\,j ^ P- 
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